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Gauge invariance and geometric phase in nonequilibrium thermodynamics 


Simone Borlenghi^ 

^Department of Physies and Astronomy, Uppsala University, Box 516, SE-75120 Uppsala, Sweden. 

We show the link between Ui lattice gauge theories and the off-equilibrium thermodynamics of a 
large class of nonlinear oscillators networks. The coupling between the oscillators plays the role of a 
gauge field, or connection, on the network. The thermodynamical forces that drive energy ffows are 
expressed in terms of the curvature of the connection, analogous to a geometric phase. The model, 
which holds both close and far from equilibrium, predicts the existence of persistent energy and 
particle currents circulating in close loops through the network. The predictions are confirmed by 
numerical simulations. Possible extension of the theory and experimental applications to nanoscale 
devices are briefly discussed. 

PACS numbers: 


I. INTRODUCTION 

Gauge theories have played a prominent role in the de¬ 
velopment of modern Physics, spanning from the descrip¬ 
tion of elementary particles to general relativity and con¬ 
densed matter theory [1-6]. However, the role of gauge 
invariance in the formulation of non-equilibrium thermo¬ 
dynamics is still poorly understood. 

The standard approach to non equilibrium transport 
processes is formulated in terms of generalised forces and 
coupled currents [7-10]. In the study of heat transport 
in nonlinear lattices [11, 12], one typically considers a 
system connected to several thermal reservoirs. 

The equations of motion of the ensemble (sys- 
tem+reservoirs) are solved numerically as a function of 
different reservoir parameters, such as temperature and 
chemical potential. 

The relevant currents are calculated in the stationary 
state and related to the thermodynamical forces by com¬ 
puting the Onsager matrix. The forces, which relate the 
dynamics to the thermodynamics of the system, are in 
general expressed as temperature and chemical poten¬ 
tial differences or gradients. The main limitation of this 
formulation consists in that it is valid close to thermal 
equilibrium. 

A more general approach was developed by Schnaken- 
berg [13]. Here a Markovian process is described by a 
master equation and represented by a graph, where the 
nodes correspond to probabilities associated to the states 
of the system and the links to transition probabilities 
(or amplitudes) between the states. Thermodynamical 
forces, that drive the probability currents in the system, 
are expressed in terms of close paths in the graph. From 
probability currents, heat, energy and particle flows can 
be obtained. The formulation holds arbitrary far from 
equilibrium and applies both to classical and quantum 
master equations [14]. 

Recently, M. Polettini has shown that the Schnaken- 
berg model can be formulated as a gauge theory [15], 
where thermodynamical forces play the role of a gauge 
potential. 

Other lines of research were focussed on the covariant 


formulation of the Fokker-Planck equation [16, 17] and on 
relating the Onsager reciprocity relations to global gauge 
symmetries [18]. 

In the present Paper, we investigate the role of gauge 
symmetries in the transport properties of a network of 
nonlinear oscillators, modelled by a generalisation of the 
discrete nonlinear Schrodinger equation (DNLS) [19, 20]. 

The recent advances on the non equilibrium DNLS [21- 
25] have shown that transport depends on the synchroni¬ 
sation and phase differences between the oscillators, shar¬ 
ing strong similarities with the well known Josephson ef¬ 
fect [26]. 

In the present work, we formalise those observations by 
showing that the coupling between the oscillators plays 
the role of a gauge field, or connection, on the network. 
The thermodynamical forces that drive energy flows are 
expressed in terms of the curvature of the connection, 
analogous to a geometric phase. This ultimately stems 
from the fluctuation-dissipation theorem and from the Ui 
gauge invariance of the DNLS. 

Our model predicts in particular the existence of per¬ 
sistent energy and particle currents circulating in close 
loops through the network, as shown in Fig. la). 

The present Paper is organised as follows. In Sec. H 
we review transport in the DNLS model. Sec. HI con¬ 
tains the main result of the paper, where we show the 
link between the off-equilibrium thermodynamics of os¬ 
cillator networks and lattice gauge theory. Here we iden¬ 
tify thermodynamical forces with lattice gauge fields and 
describe circulating currents. Sec. IV contains numeri¬ 
cal simulations that further inspect and corroborate the 
model. The conclusions, together with a summary of the 
main results and possible developments, are reported in 
Sec. V. 


II. OSCILLATOR NETWORK MODEL 

The dynamics of several physical systems, such as 
Bose-Einstein condensates, photonics waveguides, lasers, 
mechanical oscillators, spin and electronic systems [20, 
23, 27-31], can be described by the following universal 
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Figure 1: (Color online) a) Cartoon of an oscillator network 
described in terms of nodes 'ijjm and links Amn — 

The direction of the current is indicated by the black arrows. 
The phase fimn controls persistent particle and energy cur¬ 
rents along a close path P on the network (shaded grey 

area), b) Plaquette of four oscillators, described in terms of 
nodes and links coordinates. The link (m,/x) connects the 
node m to its neighbour n = m + /x in the fi direction. The 
phase /3mn, which in this notation reads [3^ (m) , plays the role 
of a gauge field, similar to a vector potential. The to¬ 
tal phase around the plaquette is the thermodynamical 
force that drives circulating currents, similar to a magnetic 
field (see text). 


oscillator model: 


Eq.(2) is the continuity equation that relates the time 
evolution of the oscillator power pm = to the ’’par¬ 

ticle” current 


imn = (4) 

between oscillators {m,n). The ensemble-averaged cur¬ 
rent, which expresses the statistical correlation between 
the oscillators, can describe various transport processes, 
such as the particle flow in Bose-Einstein condensates 
and the spin wave current in magnetic systems. 

Eq.(3) describes the dynamics of the phases 0rn(t), 
whose evolution depends on the local nonlinear frequen¬ 
cies, on the relative powers and on the phase difference 

^mn(^) — 0m (^) 0n(^)* 

In both Eqs.(2) and (3), temperature is described by 
means of the new random variable which 

has the same statistical properties as Upon ensemble 
averaging, the stochastic terms vanish and the action of 
the bath appears through 2Dm{Pm)Tm in Eq.(2). This 
source term ensures that the powers Pm are never zero at 
finite temperature. In non equilibrium steady states, Pm 
is zero in average and Eq.(2) becomes 

_ ‘^Dm{Pm)Tm + 'Yln^rnn 


^0m — ^rn{Prn)'^rn rn{Prn)'^rt 

T ^ ^ ^mn^n “k 




( 1 ) 


Eq.(l) generalises the DNLS to a network of nonlinear os¬ 
cillators with complex amplitudes 0^ = 
and the geometry specified by the coupling Amn = 

The first term on the right hand side of Eq.(l) is 
the nonlinear frequency, while Tm{Pm) = c^m^miPm) 
is the nonlinear damping, proportional to the parame¬ 
ter am- Thermal baths are described in terms of the 
complex Gaussian random variable ^m(^) with zero av¬ 
erage and variance (Cm(^)Cn(^')) = Dm{Pm)Tm^mn^{t - 
t') [22, 24, 29]. The nonlinear diffusion constant 
Dm{Pm) = ><^m{Pm)/^m{Pm) IS prescribed by the 

fluctuation-dissipation theorem [22, 29]. Tm is the bath 
temperature and A is a parameter that depends on the 
geometry of the oscillators [29]. The chemical potential 
7 ^ controls the relaxation time towards the reservoirs by 
compensating the damping. 

Physical insight is gained by writing Eq.(l) in the 
phase-amplitude representation as [22-24, 29]. 


Pm —2 [7m ^ m{Pmy\Pm T ‘^'\/Pm^^\^m (^)] 


T ‘^Dm{Pm)'^m T ^ ^ j) 
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mn 


( 2 ) 


0m — ^m{Pm) “k ^-I^[^m(^)] 

y/Pm 

^ ^ Cmn \ / COs( + Pmn)- 

^ \ Pm 


( 3 ) 


which shows that the power depends on the balance be¬ 
tween sources, losses and currents. The values of the 
currents depend on the distribution of temperature and 
chemical potential (Tm:7m)- When those are uniform, 
the system reaches thermal equilibrium, where currents 
vanish and the following equipartition relation holds: 


Pm — 


Dm{Pm)T 

^m{Pm) 7 


(6) 


The model Eq.(l) contains another current. This can 
be seen by noting that the conservative part of Eq.(l) is 
given by the functional derivative ipm = of 

the Hamiltonian [22, 29] 


^ ^ ^^ [^m{Pm)Pm “k Amn'4^m'^'n “k C.C.]. (7) 

m,n 

Computing the time evolution of PL gives the energy cur¬ 
rent [11, 21-24] 


Jmn ~ (8) 

Close to thermal equilibrium, the two coupled currents 
Eqs.(4) and (8) are related to the thermodynamical forces 
(differences of temperatures and chemical potentials) 
through the Onsager matrix, according to the standard 
formulation of non equilibrium thermodynamics [21]. In 
the stationary state, the currents are conveniently writ¬ 
ten in the phase-amplitude representation as 

Jmn ~ ‘^^mny/PmPn ^^^{^mn “k ^mn\ (9) 

Jmn ~ ^Cmny/PmPn^m {Pm ) All(^/\mn “k l^mn)- (fO) 
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From Eqs.(9) and (10) one can see that the currents can 
flow whenever the oscillators are phase locked, so that 
Amn reaches a constant value. The latter depends on 
the distribution of temperatures and chemical potentials 
trough Eqs.( 2 ) and (3). If the oscillators are not phase- 
locked, Amn fluctuates in time, so that the currents os¬ 
cillate around zero and vanish in average [23, 24]. 

Erom the fluctuation dissipation theorem follows that, 
at uniform (T, 7 ) thermal equilibrium can be attained 
only if the coupling is dissipative: Amn = ^mn[^ — 
iTm{Pm)]^ where Amn is a real matrix [ 22 ]. 

This amounts to fixing 

'’”"="(7^)^ (II) 

In this case, fdmn = Pmm is a local quantity that contains 
only the m index. Setting a nonlocal Prnm Ih^it depends 
on two oscillator indexes (or, in other terms, on the link 
between oscillators m and n) drives the system out of 
equilibrium. 

It has been recently shown by the present author that 
Pmn acts as an additional thermodynamical force that 
drives the current, in a way similar to the Josephson ef¬ 
fect [25] In particular, it permits to propagate persistent 
currents between oscillators at the same temperature, of 
from the colder to the hotter oscillator, operating the 
system as a heat pump. 

The relation between the locality/nonlocality of fdmm 
and thermodynamical forces will be thoroughly discussed 
in the next section. 


III. CONNECTION WITH LATTICE GAUGE 
THEORY 

The possibility to control transport by tuning P^nn sug¬ 
gests that it should be possible to propagate persistent 
currents = 'Ylip imrf circulating along a close path 

P in the oscillator network, as shown in Eigs.la) and 2a). 

This unusual heat circuit where energy propagates in a 
close loop, might seems a paradox, since flows, described 
phenomenologically by thermodynamical forces, are due 
to the inhomogeneous spatial distribution of some physi¬ 
cal quantities, such as concentrations and kinetic energy. 
The latter case is depicted in Eig.2b), where the currents 
generated by a temperature difference flow from the hot¬ 
ter (T+) to the colder (T_) oscillator, and there there is 
no circulating current. 

Erom Eqs.(9) and (10), it appears that persistent 
circulating currents require that the total phase 4> = 
^mneP accumulated along a close path P (or an- 
holonomy angle) is non zero. Here we shall see that fdmn 
plays the role of a Gauge field and 4> is its curvature, that 
corresponds to the thermodynamical force that drives cir¬ 
culating currents. 

Let us start from the stochastic, complex Lagrangian 
of the problem, which describes the ensemble (sys- 
tem-|-reservoirs). 
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Figure 2 : (Color online) Cartoon showing four oscillators 
(represented by the dots) coupled to thermal baths (repre¬ 
sented by the coloured squares at the corners). The direction 
of the currents is indicated by the black arrows, a) Persis¬ 
tent currents circulating between oscillators at the same 

temperature and chemical potential (T, 7 ), controlled by /3mn- 
Counterclockwise (resp. clockwise) currents are considered 
positive (resp. negative), b) Currents controlled solely by the 
differences between temperatures T+ > To > T_. Those cur¬ 
rents do not circulate, but flow from the hotter ('^s) to the 
colder (' 0 i) oscillator. 




T 'l\/ (Tm)^m(Cm'0rn ~ 


( 12 ) 


From the variation of the action 5 f Cdt = 0, one gets the 
Euler-Lagrange equations ^ -= 0, which cor- 

respond to Eq.(l). The imaginary and stochastic parts of 
the Lagrangian account respectively for dissipation and 
thermal baths. The coupling matrix Amn in general is 
not Hermitian. 

The Lagrangian Eq.(12) is unchanged by the global 
Ui rotation ipm Thus, it is left in¬ 

variant by the infinitesimal transformation Sipm ~ 
-iaipm if and only if it satisfies the equation SC = 


= 0 . 


+ St'S* + 

A straightforward calculation shows that the latter is 
precisely the continuity equation Eq.(2). In particular. 
Pm corresponds to the Noether charge and to the 
associated current. 

A crucial point here is that the presence of the stochas¬ 
tic term ^m does not break the Ui symmetry in average, 
since the statistical properties of the bath are unchanged 
by phase transformations. 

Equations and ( 12 ) is also left invariant by the more 
general, local gauge transformation 


Pmn ^ Pmn (-^n 




(13) 


Physically, this states the equivalence between the dif¬ 
ferent reference frames associated to the oscillators. In 
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other terms, the initial phases A^s of the t/^^s can be 
arbitrarily chosen, but they have to be compensated a 
by proper rescaling of the coupling to leave the dynamics 
invariant. 

All the physical observables have to be gauge invariant. 
Together with the local powers and currents {pmdrrln) 
the anholonomy angle 4> is invariant, since the sum 
— An) Vanishes along a close path P. 

The meaning of 4> becomes clear when we consider 
the simplest realisation of 2-dimensional lattice shown in 
Fig. lb), consisting of four coupled oscillators and called 
plaquette. More complicate lattices can be constructed 
by joining several plaquettes. By adopting the notation 
commonly used in lattice gauge theories [1, 2], we de¬ 
scribe the plaquette in terms of its nodes m, and links 
between the nodes along the coordinate directions (/i, u). 

We label with the indexes (m, /i, u) the directed link 
between oscillator m and its neighbour along the p or u 
direction. Using this notation, the coupling matrix reads 
where P^{m) corresponds to the phase 
Prnn th^t connccts oscillators m and n = m-\-p. By intro¬ 
ducing the finite difference gradient operator V^/d(m) = 
and upon defining /3-^{mPp) = 

[2], the anholonomy angle around the plaquette reads 

^ - V^/3^(m), (14) 

which is the discretised version of the curl operator. Us¬ 
ing the same notation, the gauge transformation Eq.(13) 
reads ^^(m) ^ ^^(m) - V^A(m). 

Here /3^{m) plays the role of a gauge field (connection) 
on the oscillator lattice, analogous to the vector poten¬ 
tial of electromagnetism. Similarly, the total phase 
is the curvature of the connection, analogous to the Fara¬ 
day field tensor. This is equivalent to a geometric phase, 
that depends only on the geometric properties of the path 
and not on the dynamics. On the other hand, the phase- 
differences between the oscillators in Eq.(2) can be 
seen as dynamical phases, that depend on all the pa¬ 
rameters (T^, 7 ^,/3^(m)), thus on the system dynamics. 
They vanish on close loops and therefore do not con¬ 
tribute to current circulation. 

At thermal equilibrium, the phases are local quanti¬ 
ties Prn that obey Eq.(ll). This guarantees that the an¬ 
holonomy angle 4> is zero, since Eq.(14) reduces to the 
sum — Pn) which vanishes on close paths. On 

the contrary, one can have a situation where there the 
system is out of equilibrium, but 4> = 0. In this case 
the local currents j^n between neighbouring oscillators 
do not vanish, but their sum around the path is zero, so 
that there is no net circulating current. 

Gauge invariants can be in general expressed in terms 
of Wilson loops, i.e. path-oriented exponentials W = 
Pexp{J]]^^p /3^}, proportional to path ordered products 
of A^. The latter can be also identified as the the gauge 
field instead of The Noether currents are gauge in¬ 
variants given by the coupling of the gauge field with 



Figure 3: (Color online). Time evolution of the system at zero 
temperature and phase ydmn- The off-equilibrium dynamics is 
given by a nonuniform distribution of chemical potential /im, 
m = 1,..., 4. Panel (a) shows the dynamics of the powers pn, 
that reach constant values in the stationary state. The inset 
shows a detail of the dynamics in the stationary state, where 
one can see the difference between the powers. In panel (b) 
the time evolution of the particle currents E displayed. 
The currents are nonzero, but their sum along the path 
(displayed in orange tones) vanishes as expected 


the oscillator wavefunctions (V^m^V^n), that play the role 
of ” matter fields” [2]. 


IV. NUMERICAL SIMULATIONS 

To corroborate the model, we turn now to numerical 
simulations. We consider the off-equilibrium dynamics 
of a plaquette, made of the coupled oscillators ^ = 
1,..., 4 shown in Fig. 2a). 

The oscillators have the same nonlinear frequencies 
and damping rates, respectively uJn{Pn) = x (1 + 2 pn) 
and Tn{pn) = <^^n{Pn) (l^ uuits where ks = 1). Each 
oscillator is connected to an independent Langevin bath. 
The baths have temperatures and chemical potentials 
(T^, 7 ^), and the same coupling strength D = 2a. 

Throughout all this section, for simplicity we adopt the 
notation of Sec.II and Eig. 2a), with oscillator indexes 
(m, n). The coupling between oscillators (m,n) reads 
thenA^^ = (7(1 — ia)e^^'^'^. The anholonomy angle, or 
gauge field curvature, reads = ^12 + ^23 + /d 34 -h Pai. 

The equations of motion were solved numerically using 
a fourth order Runge-Kutta algorithm with a time step 
dt = 10“^. The other parameters adopted in all our 
simulations are = 1 and C = 0.1. All the parameters 
are expressed in model units. 

We consider first the case at zero temperature, and 
relatively high damping a = 0.1. Eor simplicity only the 
particle currents are reported, since the energy currents 
have the same profiles up to a scaling factor. 

Eigs.3 (a) and (b) show respectively the time evolution 
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Figure 4: (Color online). Time evolution of the system at zero 
temperature and chemical potential. The system is driven 
off-equilibrium solely by the phases /3mn- Panel (a) shows 
the dynamics of the powers pn, that reach constant values in 
the stationary state. Panel (b) shows the time evolution of 
the particle currents • Their sum along the path (dis¬ 
played in orange tones) does not vanish, since the anholonomy 
angle 4> is nonzero. 


of the powers Pm and of the particle currents (j^^, J^). 
Here (3^^ = 0 and the following values of chemical poten¬ 
tial where chosen: (71 = 0 . 13,72 = 0.15,73 = 0 . 18,74 = 
0.11). After a transient regime the system reaches a non 
equilibrium stationary state where powers are different 
and local currents are constant in time. However, as 
expected from our model, the circulating current along 
the plaquette (denoted in orange tones) is zero. 

Figs. 4 (a) and (b) again display respectively the time 
evolution of powers and particle currents. In this case 
all chemical potentials are set to zero, while [3mn has 
the following values: {P 12 = 0 . 257 r ,;523 = 0 . 577 r ,;534 = 
O.SStt, = 0.327r) These phases drive the system in a 
non equilibrium state where all the currents, including 
the sum around the plaquette are nonzero. 

Next, we consider the more general case where both 
temperature and chemical potentials are nonzero. In the 
following simulations the observables were time averaged 
in the stationary state over an interval of 3.5 x 10^ time 
steps and then ensemble-averaged over 100 samples with 
different realisations of the thermal field. 

We have taken a smaller damping parameter a = 0 . 02 , 
and we have set (/ 3 i 2 = 0 . 57 r ,/323 = 1.27r,;534 = O.Stt). 
The curvature reads then ^ = 0 . 497 r+/ 54 i. The dynamics 
was computed as a function of / 54 i, in the presence of 
uniform temperature and chemical potential (T, 7 ). 

Fig. 5 (a) and (b) shows respectively the powers and 
currents vs / 34 i. The computation where performed with 
bath parameters T = 0.3 and 7 = 0 . 021 . One can see 
clearly that persistent currents circulate through the sys¬ 
tem. When the currents change, energy is redistributed 
through the system and also the local powers change. 
This testifies that the current circulation corresponds to 



Figure 5: (Color online) a) Oscillator powers and b) total 
currents vs / 34 i, computed using the bath parameters T = 0.3 
and 7 = 0.021 model units. The current is magnihed by a 
factor 3 for better visibility. The solid lines are guides to the 
eye. 



Figure 6 : (Color online) Oscillator powers (a) and currents 
(b) vs chemical potential 7 , computed for ^ 4.1 = 0 . 37 r and 
T = 0.3. Powers and currents grow with 7 , as the energy of 
the system increases. The current is magnified by a factor 
6 for better visibility. The dashed black lines are hts with a 
Fade function (see text). Since the powers have very close 
values, for better visibility in panel (a) only the ht for pi is 
displayed. 


a real transport process through the system. Apart from 
their different magnitudes, the currents have simi¬ 

lar profiles. 

By increasing (T, 7 ), the powers and energy of the sys¬ 
tem grow. One expect that this leads also to an in¬ 
crease in the currents. This is indeed the case when 7 
increases and T is kept fixed, as one can see in Fig. 4 (a) 
and (b). Here powers and currents vs 7 are displayed, 
with parameters T = 0.03 and / 34 i = 0 . 37 r. Both powers 
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Figure 7: (Color online) Oscillator powers (a) and currents 
(b) vs temperature T, computed for / 34 i = O.Stt and 7 = 
0.021. Powers and currents decrease with temperature T, 
since thermal fluctuations reduce the coherence between the 
oscillators. The current is magnified by a factor 6 for 
better visibility. The lines are guides to the eye. 


and currents are fitted with a Fade function ^ 9 ( 7 ) (resp. 
(ao + ai7 + a27^)/(l - «37)- 
However, when temperature increases and 7 = 0.021 
remains fixed, the situation is reversed: powers and cur¬ 
rents drop fast and reach an asymptotic value around 
T = 0.3. This is most likely due to the fact that temper¬ 
ature does not increase only energy, but also phase fluctu¬ 
ations. Those break the exact synchronisation and tend 
to suppress transport processes. This feature shows the 
importance of coherence in energy transport processes of 
DNLS systems. 

V. CONCLUSIONS 

We have demonstrated that the off-equilibrium ther¬ 
modynamics of a large class of oscillating systems can be 
described by a Ui lattice gauge theory. The gauge fields 
correspond to the complex coupling between the oscil¬ 
lators. The thermodynamical forces that generate cir¬ 


culating currents are the curvature of the fields. Those 
forces are thus geometrical properties of the ensemble 
(system-breservoirs). 

The main prediction of the model, the circulation of 
persistent currents along close loops, is confirmed by nu¬ 
merical simulations. This result is robust within a large 
range of temperature of physical parameters and hold 
far from thermal equilibrium. We find that transport in¬ 
creases with chemical potential but is suppressed by ther¬ 
mal fluctuations that tend to destroy the phase coherence 
of the system, a behaviour coherent with previous studies 
on the DNLS. 

This work open the path to several possible develop¬ 
ments. In particular, the off equilibrium thermodynam¬ 
ics of system with other kind of gauge symmetries [ 2 ] 
could be studied in a similar way. The transition be¬ 
tween coherent and incoherent transport, where currents 
are carried mainly by the modulation of the powers and 
the phase are completely incoherent, should also be in¬ 
vestigated. 

The role of close paths in our formulation suggests a 
link with Schnakenberger theory, and further investiga¬ 
tion is needed in this direction. 

Concerning possible experiments, current circulation 
could possibly be observed in spin system with the 
Dzyaloshinskii-Moriya Interaction [32], characterised by 
an asymmetric propagation of spin waves. The dynam¬ 
ics of those systems can be phenomenologically described 
by means of Eq.(l) [29] with a complex asymmetric cou¬ 
pling term 7 ^ ^nm- This could allow for a nonzero 
anholonomy angle on close paths. 
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